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The Incidence of Light upon a Transparent Sphere of Dimensions 

Comparable with the Wave-Length. 

By Loed Bayleigh, O.M., F.E.S. 

(Received April 4, — Read April 21, 1910.) 

In a paper, now nearly thirty years old,* I applied Maxwell's equations 
of the electro-magnetic field to investigate the disturbance produced by an 
obstacle upon plane waves of light which travel through a medium other- 
wise uniform, giving particular attention to the case where the properties 
of the obstacle differ but little from those of its surroundings. The difference 
may consist in a variation of K~ the specific inductive capacity, or of 
yu,— the magnetic capacity, or of both; but it was shown that the last 
supposition leads to results inconsistent with observation, and that the 
evidence favours the view that jjl is to be treated as invariable. Denoting 
electric displacements by /, g, h, the primary wave was taken to be 

h = e int e ikx , (23) 

so that the direction of propagation is along x (negatively), and that of 
vibration parallel to z. A/x being omitted, the electric displacements 
(fh 9i, h{) in the scattered wave, so far as they depend upon the first power 
of AK, have at a great distance the values 

fl , gi , h = J ^(q, ^ - ^) , (35, 37, 38) 

4tirr \r 2 r 2 r 2 I 

Ao AK-V" ttr dx dy dz. (36) 



in which P = 



In these equations r denotes the distance between the point (a, /3, 7), 
where the disturbance is required to be estimated, and the element of 
volume of the obstacle dx dy dz. 

It is evidently implied that in a direction (a = 0, /3 = 0) parallel to that 
of primary vibration, the scattered vibration vanishes, whatever may be the 
size and shape of the obstacle, so that any light which may appear in this 
direction must depend upon powers of AK higher than the first. 

If AK be uniform throughout the obstacle, it may be removed from under 
the integral sign in (36), and if the value of K for the external medium be 
taken as unity, we may write 

KAK" 1 = -(K-l). 

* 'Phil. Mag., 5 vol. 12, p. 81, 1881 ; 'Scientific Papers, 5 vol. 1, p. 518. In the second 
term of Equation (32) and following, AK™ 1 should read Afi~\ 
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If the obstacle is very small, h and r may be treated as constants in the 

integration, and 

p _ -(K-l)h e- ikr .T, (40) 

if T be the whole volume of the small obstacle. Using this in (35), etc., 
we get 



TrT(K-l) , 



\ 2 r 



ocy f3y a 2 -f /3 



2" 



/i, 9x, h = " M \7 > * (nt - kr) ~H> -H> ^r~ , (41, 42, 43) 



rt->^ ,™J ^>^ 



A, being the wave-length (27r/k). These equations apply, whatever may be 
the shape of the small obstacle ; and, of course, they undergo no change if 
the shape be spherical. In that case it was shown that the restriction upon 
the value of K— -1 may be removed, provided we replace K— -1 by 

3(K-l)/(K + 2). 

We have seen that the scattered light vanishes to a first approximation 
in (K— 1) for the direction parallel to z. In the case of the sphere it was 
shown that for this direction we have to a second approximation for a small 
obstacle 

/„ = 7r JyK~-I) e i( n t-kr) m * ^ /go) 

with g 2 = 0, h 2 = 0. In (60) K denotes the radius of the sphere, and r is 
measured from the centre. " Comparing (60) and (41) we see that the 
amplitude of the light scattered along z is not only of higher order in AK, 
but is also of the order h 2 W in comparison with that scattered in other 
directions. The incident light being white, the intensity of the component 
colours scattered along z varies as the inverse eighth power of the wave- 
length, so that the resultant light is a rich blue." It is obvious from (41) 
that in general the intensity of the scattered light varies as X" 4 . 

" There is another point of importance to be noticed. Although when the 
terms of the second order are included, the scattered light does not vanish 
along the axis of z, the peculiarity is not lost but merely transferred to 
another direction. Putting together the terms of the first and second orders 
we see that the scattered light will vanish in a direction in the plane of xz, 
inclined to z (towards + x) at a small angle 0, such that 

A"K" 7AR 2 fr 2 "R 2 

^ = ir^l" = (K ~ 1) ^-" (61) 

Some experiments in illustration were then described. 

Eeturning to the case where (K — l) 2 is neglected, we consider a spherical 

obstacle of any radius E, whose centre is at the origin of co-ordinates. 

From (23), (36) 

P = — (K — 1) e int \l\ e, ik (x-r) dx dy dz, 
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r being the distance between dx dy dz, and the (distant) point at which 
fhffh^i are to be estimated. It is evident that, so far as the secondary ray 
is concerned, P depends only upon the angle (^) which this ray makes with 
the primary ray. We will suppose that % = in the direction backwards 
along the primary ray, and that % = it along the primary ray continued. 
The integral may then be found in the form 

fccosiv J Jl ( 2 * K cos ix cos 0) cos2 # ; ^ ' 

Tq being the distance of the point of observation from the centre of the 
sphere. If we expand the Bessel's function and drop the suffix in r as no 
longer required, we get 



3 I 2.5 2.4.5. 



m 6 , m 8 "^ 



2.4.6.5.7.9 2.4.6.8.5.7.9.11 



7 

(49) 



in which m is written for 2&Ecos |-^. From this fi,gi, hi follow at once by 
(35), etc. Along the continuation of the primary ray cos |% = 0, and P 
reduces to 



4 7T 



R 8 (K-l)e*<»'-* r >, 



3 

as was to be expected. It is to be observed that in this solution there is no 
limitation upon the value of K if (K— l) 2 is neglected absolutely. 

Prof. Love, in a valuable paper on the " Scattering of Electric Waves by 
a Dielectric Sphere,"* has treated this problem by a different method, 
limited on the one hand bo the sphere, but on the other applicable whatever 
may be the value of (K — 1), and (so far as the general analytical expressions 
are concerned) whatever may be the size of the sphere. In the 1881 paper 
I had treated in this way the problem of a dielectric cylinder. From these 
expressions Love deduces first and second approximations applicable when 
the radius of the sphere may be treated as small in comparison with the 
wave-length. It will be of interest to compare these approximations with 
those already referred to. 

In Prof. Love's notation the electric forces are represented by X, Y, Z ; 
and, so far as relates to the medium outside the sphere, they may be 
identified with /, g, h used above. The co-ordinates of the point of observa- 
tion are denoted by x, y 9 % instead of by #, /3, 7. Again, Love supposes the 
direction of propagation to be along — z (instead of —x), and the direction of 

* ' Lond. Math. Soc. Proc.,' vol. 30, p. 308, 1899. 



28 Lord Rayleigh. Incidence of Light on a Transparent [Apr. 4, 

vibration to be parallel to y (instead of z), changes represented by taking one 
step backwards in the cycle x, y, z } x. Thus (with omission of a constant 
multiplier which runs through) the primary wave is represented by 

Y = e *{ct+z) m (X) 

Corresponding to this, the complete expressions, so far as terms of order MR 5 , 
for the forces in the scattered wave at a great distance are found to be 

yz> 



(X v 7 \ „ K-17£ 2 RV'* (c *" r+R) / xy x 2 + z 2 



rp2i 



I \18 5K + 2/J 



+(K-i)iiLi- — Vo,--, y 

30 r \ r r 

2K + 3 Qr \ ' r' r 

+ 9T . , Q Q ^ -, '', ""I. (XLH) 

ZJ1 + 3 or r J \r r rj 

So far as the first approximation involving WW, this agrees mutatis 
mutandis with (41, 42, 43), since T = |-7rPt 3 and h = 27r/X. 

We may next consider the application of (XLII) to the secondary ray in 
the direction of primary vibration (x = 0, z = 0, y = r) } for which the first 
approximation vanishes. We have 

(X, Y, Z) = ^ ^ (0, 0, 1) 4- ^-^ (0, 0,-1) 



30 r 



2K + 3 



_ (K-l) 2 FE 5 ^^- y ) /n ft n 
"" 15(2K + 3)r W u '^' 



and this agrees with (60) when we introduce the supposition that (K— 1) is 
small, so that 2K + 3 in the denomination may be identified with 5. 

I have already commented upon this agreement ;* but recently in pursuing 
the comparison I came upon a discrepancy. It is evident from the theory 
already given that a great general simplification should attend the supposition 
that (K — 1) is very small. Introducing this into (XLII) and omitting for 
the moment the factor 

pikifit— r+H) 

(K-1U 2 K 3 - , 

r 

we find for the three last sets of terms 

k 2 R 2 z 



xy x 2 -\-z 2 __yz\ 

/via /ltd ty>A i 



* 'Scientific Papers/ vol. 1, p. 536, 1899. 



1910.] Sphere of Dimensions Comparable with Wave-Length. 29 



showing within the bracket the same dependence on x, y, z as does the first 
set. Hence altogether, expanding the exponential e ikR , we may write in this 
case 

(A, IjZj) — ^ I r 2 » -2 » ..2 



r 



.2 



x[l + P^(i-n+f) 



PE 2 ^ 



o r J 



This should correspond with (35), (49) of the former investigation if in the 
latter we stop at the term in m 2 where 



z 



m 2 = 4 WS? cos 2 1- % = 2 PE 2 1 + - ) , 



so that 



P= - 



4ttE 3 



(K-ij^r. 

3 I 



rj 

5 \ r i 



It will be seen that in all respects except one the agreement is complete. 
For example, i has disappeared otherwise than as contained in e^ ct ~ r \ But 
there remains a numerical discrepancy, since (2— rf + "§) * s no ^ ec L ua l to — -J. 

For some time I was at a loss to account for this discrepancy, until I 

noticed that the omission of a term -^ WS? in Love's equation (XXXIV) 

would produce harmony, and this led me to scrutinise equation (XVIII), 

where (I believe) the error lies. It appears to me that these equations 

should be 

jjceikct r t?2 



4>n = 



n (n+l) yfr n (kR)L.2n + 3 



\: 



E2 ot»+l(^) rfV " +1 



dx 



,2n+l 



2n-l 



^r„-l(*KW 



d /V„-i' 



dx\r 2n 1 /_ 



X* 



ik e lkct 



n(n + l)yfr n (k'R) 



W 



.2n + 3 



*-C*B)^ 



Kxvin*) 



ryi. 



2re+l 



2t&-1 



^n-l(AR) 



d /Vn-l\ 

dy\r 2n ~ 1 i 



the change consisting in the substitution of (2w + 3) for (2^ + 1) in the 
denominators of the first terms of <f) n , % n . We should then obtain 

^ = __ 1 ^ gifatf ^ simply, 

instead of (XXXIV), which involves also the factor (1+ T 2 ¥ & 2 B 2 ). This is a 
matter to which I shall presently return. 

In view of numerical calculations, it may be well to complete the statement 
of Prof. Love's results, and to show that, especially when corrected as above, 
they admit of considerable simplification. 
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Tor the functions which occur in (XVIII) we have Y = 1, and for n 
greater than unity 

V — (iJcr) n p I z\ (YIY^ 

n ~1.3.5...(2tt-l) in W' {J ^ iy) 

P w denoting as usual Legendre's function, so that V n is a solid harmonic of 
degree n. For ^r n we have 

\V M) I V 

— 1 ^ 2 , rf 7f 

"" 2.2ti + 3 2.4.2^ + 3. 2ti + 5 2. 4. 6(2^+3) (2 ti + 5) (2 ti + 7) 

+ ..., (IV)* 
and the functions ^ w (77) satisfy certain sequence equations which may be 
written 

J^fnO?) = - 2 -^3*.+iW = (2»+l){^-ifo)-*„fo)}. (VI) 

We have also to deal with another function, E^ (77), such that 

E n (rj) = (-j»1.3.5...(2n+l)(±4)*—> (VII) 

and the functions E w (77) satisfy the same sequence equations as the functions 
-^r n (rj). It will be seen that —i^r n (v) is the imaginary part of E n (r)). As 
included in (VII), 

E (77) = — , Ei (n) = 3 *±^ erK (XXXVI) 

77 77^ 

This being understood, the disturbance in the secondary wave, so far as it 
depends upon <j> m % n , may be written 



X n = M J JE»(Ar)(y^-^ 



1ST. 



+ -T 



-( w + l)K. 1 (Ar)^ + — ^^^iWil 



with two similar equations for Y n , Z n . The complete expressions for X, Y, Z 
involve a summation with respect to n which takes the values 1, 2, 3, etc. 

All is now denned except the constants M n , N" n , which depend upon the 
boundary conditions to be satisfied at the surface of the sphere. For the 
present case of a dielectric sphere of refractive index ^K, or Jtf/Jc, 
Prof. Love finds the equivalent of 

^-i(fcK)-*P^>*.(fcK) 
M M = y» { -7 , (XXIII) 

-E.-X (*K) + ^rlj-ffi E re (fcK) 

Tn (« K) 

* This series is the same as that which occurs in the expression for Jb+j(ij). 
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-KE,-! QeR) + {(K - 1) ^— + *f~}^\ E„ (*E) 

L 2^ + 1 i|r n (iE) J 

(XXIII Us) 

which completes the analytical solution. 

Eeturning now to (XV1IP) we see that the first terms dV n+1 /dx, 
dV n+ i/dy are solid harmonics of degree n. Also Yn-i/r 2 * 1 " 1 is a solid 
harmonic of degree — n. When differentiated with respect to x or y, it 
becomes a solid harmonic of degree —n--l, and when further multiplied by 
r 2w+1 , a solid harmonic of degree n. Hence both terms in (XVIII*) are solid 
harmonics of degree n, and it will appear (as might, indeed, be anticipated) 
that they are of the same form. Writing p for z/r, we have from (XIV) 

~dx~ - 1.3. 5... (271+1) j^ 1 )^^)-^ -^ j- • (*) 

Again 

far**-* 1.3.5...(2»-3)1 B ^»-iW- ^-^T/- W 

Now (Todhunter's 'Laplace's Functions,' pp. '43, 46) the quantities 
within { } in (a) and (5) are equal to one another and to —dP n /d/x,, so that 

,.2»+i d_ Yg-i _ _ (2^-l)(2n + l)^V M+1 

J R 2 i /jw^n+l r 2M+1 , /7T>X d V»_i 

2n + 6 dx 2n — l T v 7 (fer 2n+1 

- "& F~| 2^+1.2^ + 3 ^ l(/;R) + ^- l(Z;R) 

_ dY n +i 2 n+1 , ...p. 

— ^" Hyr" *"<**>' 

by (VI). We thus obtain the greatly simplified forms 



, _ (2w+l)^-V« f ^V w+1 _rfV B+1 - 



(d) 



in which E does not appear at all. 

Eeferring to (XIV) and introducing the value of P n (s/r), we find 

V = 1, V 2 == -£ #» (2 a»-aj8-y8), -^ 

V x = fe, V 3 = -g!{2a8-3« (^ 2 + 2/ 2 )}.J (XXVT) 

Thus in the case of n = 1 we have 

^ = ~-\%kyz ikct y 
without any approximation depending on the smallness of E. 
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For further reductions it will be convenient to introduce a modified form 
V n , viz, l) m denned by 



Y, = 



(iky it 



(ikr)* 



n 



1.3. 5. ..(2*1-1) 1.3. 5. ..(2»-l) n \W' K } 



We have now to consider the solution expressed in (XV), and we may at 

once introduce the simplifying condition that r is very great. In this case 

f 
rom (V±l) 



Tfi, n (kr) = i n ,1 .3.5 ... (2^ + 1) 



e 



'-ikr 



{hr) n 



+1 > 



(XXVIII) 



and we find from (XV), (d) 



.A 



n. 



n . 7i+l . hr n * 1 



M n (y~-2 



am *<£ wy«M. i * y 



clz clyj dx 



■K 



w 



r- 



dx r } 



Tfbl 



dy __ 



(«) 



and in the cyclical changes by which Y and Z are deduced from X, dXJ n +ifd&, 
dJJ n +i/df are to stand unaltered. 

For systematic calculation it is best to express U w+ i in terms of P. 

As in (a) 



dV 






dx 



so that we may write 



»+i - % T to-i{( n + l)'P n + l —p'P' n+1 } = —a?r n ~ 1 P , n ; 



(/) 



dU 



w+l 



<ZU 



?» + l 



(XX 



<% 






(3) 



The differential coefficients which appear in (e) are now readily formed. 

e find 






AW 



w 



Y 



(• 



» 



)» + H2rc + l)e* (c *~ y) 
n.n+1 .Jcr 



M 



» 



*,P' _L P" 

- mx n -| 5 x w 



» 



AW 



// 



w 



(*) 



Z 



B 



(-)»-"(2*+l)<**< rt " ) g [M n P / B + N n {-MP / - + (l-A* s )P"»}]. 0) 

-«, . w 4- 1 . kr r 



The { } in (j) may alternatively be written fxY n —n(n + l) P„ ? in virtue of 
the differential equation satisfied by P w . 

It should be observed that the above values satisfy xX n +yY n +z Z n = 0, as 
of course they ought to do. 
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At this stage it may be convenient to write down the values of P'» for tho 
earlier values of 'n. 





P' 


n — 1 

n = 2 

n = 3 

^ = 4 


1 

3yU 

0,0.0 o 

2 2 

3.5 , 5.7 3 

2 ^ 2 * 



Thus from (h), (i), (/), 



Xi, Yi, Zi 



Q gik(ct-r) 



2 . fo- 



^2; JL 2, <^2 = — 



5 e ik(ct-r) 

2.kr 



'-1 Ni , -|*Mi + N x ^1 -£j , | (Mi -/.Ni) 



2{^M a +N-,(l-2/t a )} 



A9 = 



T» = 



1 e iHct-r) xy r /l 3.5 , 



*7 J.k{ct—r) 

4.kr _ 



MM 3 -ji-^Hsfl-^ 



2 



+ N 



y 2 \/3 

,1,2 



.5 „ 1\ 



2 M 2 



- 6 "'J- 



} 



z» = 



'3 



.X.4 = 



Y 4 



^ e ik(ct-r) y 

i.kr * r 

^ 9 e ik(ct~r) X y 

4 . kr r 2 



5 o 1\ . ™- ill 15 2 



^'i^-^/^^u"^^' 



9 



M4f l~¥ / " 2)+N4(3/ "~ l4/ " 3) 



9^(rf-r)r" j 7 3 32/3 3.7 a \ 

-4^-L M4 1-2^ + ^ -?\2-irn. 



2 



z -- 9 -^NI^)W-",< + 



3.9 



3" 



T* % 



Our expressions naturally simplify when we limit ourselves to secondary 
rays which lie in the principal planes. No real loss of generality is thereby 
incurred, since whatever may be the direction of primary vibration, that 
vibration may be resolved into two, respectively in and perpendicular to the 
plane which contains the primary and secondary rays. Thus, if y = 0, that 
VOL lxxxiv. — A. 1) 
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is for secondary rays whose direction is perpendicular to that of primary 
vibration, we have generally X n = 0, Z n = 0, and 

n.n+l.kr l 

Again, if a? = 0, that is for secondary rays lying in the plane which 
contains the directions of primary vibration and propagation, X n = 0, and 

n . n + 1 . fcr 

while Z w retains the form (/). The direction of vibration is in the plane 
x = 0. and is, of course, perpendicular to the secondary ray. Its magnitude 
is given by 

y^-zY^ = (-n2»+1)^ [Mn P' re + N. {/.P'. -» (» + 1) P.}] • (0 

r ?i.. 71 + 1 . At 

It will be seen that (k) and (Z) differ only by an interchange of M», ST„. 

It may be remarked that when p = + 1, (A;) and (7) should become 
identical, except as to sign. If /m = — 1, the signs are the same ; 
if ^ = +1, they are opposite. Thus when p = +1, the [ ] in (J) becomes 
;P' n (M^HrlsTi) simply. And in this case by the differential equation, 

"When n is even, P n = 1, and then P' n = .±|-n (?i+l) ; when w is odd, 
P„ = ±1, and then P n ' = ^?i(7i+l). 

Accordingly, (I) gives in the case of fi — ±1, 

r At 

qp(|N 1 + fM 2 + lJN T 3 + fM 4 + ...E (™), 

the upper sign in the ambiguities relating to fi = +1. 

Arranged in powers of /x, the general expression from (Z) as far as M 5> 1ST 5 
inclusive is 
ffZ-sY _g^g-^) 



r 



A;?* 



-M1--N2--M3 + — N 4 + ^M 5 

_2 2 8 8 16 






(*). 
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The corresponding expression for Y from (k) is derivable by interchange of 
M and N. 

When fx = Q, we have in (I) to consider 

In this when n is odd, V n (0) = and P' n (0) takes in succession the values 

' ~2' 271' ~2.4.6' 2.4.6.8' ' 
while when n is even, P'„ (0) = 0, and P„ (0) takes the values 

p-j-i X X . o X , o . i.o.O./ j 

L J ' ~2' 274' ~ 27*76' 2717678' 
Using these in (I), we find for /jl = 



r kr 



l Ml -f N2 -2-4 M3 + 0^ + M M5 



X t) . <J . -x-r X D , O . -»*• , X / . O • t) . / vr , 

. j\j q — ,_ — — ^ — jyi^, _j_ __ — — ^ — __ JM 3 -f- . . . 



2.4.6 2.4.6.8 ' 2.4.6.8 



(*) 



•and the corresponding value of Y in (k) is obtained by interchange of M 

and N". 

We now pass on to consider the values of M n , N w as given in (XXIII). 
By (VII), when 77 is small, "E n (v) is of the order rj~^ n+l \ so that l$ n is of the 
order rj 2n+1 . M n is two orders higher, viz., 7] 2n+3 , since the numerator is 
itself a small quantity of the order rf. When 77 is small, the most important 
terms in X, Y, Z depend upon TS\ and are of order rf or E 3 . If we are 
satisfied with an approximation as far as E 5 inclusive, we have only to 
consider Ni, Mi, and N 2 . For E 7 we need Ni, M^ N" 2 , M 2> N" 3 , and so on. 

Let us begin by supposing that (K — 1) is small, and retain only the first 
power of this quantity. For example, in *fy n (A'R), we may write 

<f n (jfc'R) = ty n (JcR) + yjr'n (AR) . (*' - k) E. 

We find, omitting for brevity the argument in ty n (k~R), etc., 

^"n-1 — «„ T -1 fJfn + 1 &E (^n-l^'n — f^Vl) 
N n = (K- 1) \ tl±l L , (q) 

The denominator in (p), (q) appears at first sight complex, since E is so 

D 2 
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If, however, we separate real and imaginary parts, writing in accordance with 
the definitions (IV), (VII) 

E B (rj) = ¥» ( V ) - % f n ( v ), (r) 

we see that 

from which the imaginary part has disappeared. In this case accordingly 
M w , IsF^ are wholly real. 

The denominator (s) admits of simple evaluation. By (VI) 

Now, if Y n be a solid harmonic, the same in both cases, Y n ^r n and Y n E n 
both satisfy the equation (y 2 -f k 2 ) = 0, and hence by Green's theorem 
(' Theory of Sound/ vol. 2, p. 252) as applied to the space included between 
two concentric spheres, we find that, 

where C n is independent of rj. To determine C^ we may suppose tj very 
small, when 

fn = i, f n = __2_, e. = (2^i) i8 - 3 -^r i)3 . 

Thus C„ = l 2 .3 2 .5 2 ...(2ra+l) 2 , 

and f»-i E »-^ E »-i = l 2 .3 3 .5 3 ...(2w-l) 2 x(2»+l)^- 2 "- 1 . (0 

It will be understood that (t) is true without approximation as regards r), 
if (K — -l) 2 is neglected. 

It is now a simple matter to calculate M n , N" n to a moderately high power 
of 7). So far as rf inclusive, we have 

1 9 V 5 10.5.7/' 45 V 7/' 

2 5.15 V 7/' 35.45 * 

7 ^?; ir 

/ » oo . *±u 

Til 
nus 

e ik(ct-r) y rg g 

Z = Zi -f- Zg -f- Zg '= — — . - -™ M.3 *—* — xv 2 



It r 



.^5 £ 



9 
O 






1910.] Sphere of Dimensions Comparable with Wave- Length. 37 

in which f Mi— -§-^"2"= 0, so that //, becomes a factor of the whole expression 
as we know it ought' to do. Also 



X 



e ik(ct~r) X y 



3 Ni + ^M 2 h4n3+m(5N 2 ~^M 3 )~^ 



hr rL2 i 2* , 8 or A 4 0/ 2.4 



In X the multiplier of xy should be the same as the multiplier of yz in Z. 
In order that this may be so, we must have M2 = \ N3, a relation satisfied by 
the tabular values. In the expression for Y, the part which does not 
contain the factor l—y 2 /r 2 y or (x 2 + z 2 ) fr 2 , vanishes in virtue of the relations 
already mentioned. Finally, arranged in powers of 77, or &R, 

y y 7 - ^ (rf - r) ^R 8 (K--i) / vy % 2 +z 2 __y^ 

in agreement with (35), (49) so far as the present solution extends. The 
former method is incomparably the more appropriate in this particular case, 
giving the result for the sphere to any degree of approximation in k~R } as well 
as admitting of application to obstacles of other forms. The general con- 
clusion that (0) vanishes when (K — l) 2 is neglected would probably be difficult 
to establish by the present method. 

Abandoning the restriction as to the smallness of (K — l), we find from 
(XXIII) as far as E 5 inclusive 



1 3(K + 2) t. * 



6 KFB 2 ~\ 
5(K + 2)J ' 



Ml = , N 2 = ___2__; ( M ) 



where it will be seen that to this order M, N are still real. The corresponding 
values of X, Y, Z are given by 

l.kr \ <r r 2, r 2 j 2 .kr \ r r 

2.kr L t s r \ r 2 / r \ r 2 ) J v J 

This agrees with Prof. Love's solution, except in the first set of terms, as 
already mentioned. 

The special case where the secondary ray lies in one or other principal 
plane is given by (n) and its derivative, when the approximate values of M, N" 
are introduced. 
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When &R, or tj, can no longer be treated as small, we must calculate M 
and N" from the general formulae (XXIII) and (XXIII bis). Separating E 
into its real and imaginary parts, we see that the denominator of M» is 
equal to 

- ¥»-i (v) + %^? Vn (v) + i x numerator. 

Yn(v) 

Also that the numerator for ~N n exceeds the numerator for M n by 

(K^l)yfr n - 1 (7 ] )-(K-l)j^^f n (7 ] ); 

and that the real part of the denominator of N n exceeds the real part of the 
denominator of M n by 

(K-lJ-^t.W-lK-l)?,..^ 

J!i 71 -j" J- 

while the imaginary part of the denominator of ~N n is i x the numerator of 
N n . When the numerators and denominators have been computed, the values 
of M n> N"„ must be reduced to the standard complex form. 

In selecting special cases for calculation we notice that 77 and rj' (JcR and 
&'K) are both arbitrary. With the view of reaching something distinctive, it 
seemed desirable to choose a refractive index Qc jlc) not too small. For this 
purpose V \h = 1*5 seemed suitable, as large enough, and yet not outside the 
experimental range. This corresponds to K= 2*25, By further choosing 
97 = 1, 7] = 3/2, 7) = 9/4 we make the tabulation of ^ serve a double 
purpose. 

The most obvious procedure for calculating E w (??), equivalent to 
^n (??) — ityn (v)> i g ^° use the sequence formula (VI), 

£ 71- -}- I . 2 71 -\- o ( ~rj xp % /^7T\ 

w +i — 5 \Ji«» — ^-ij- ? V vl / 

7) 

starting from E and E, as given by (XXXVI). So far as regards the real 
part, represented by "^P, this course is satisfactory; but the same cannot be 
said of the application to i/r. As n increases, the successive i/r's tend to 
equality, and any errors are multiplied at the next step by the large factor 
(2n+l) (2n+3). So much is this the case that although we use seven-figure 
logarithms throughout, even the first figure in i/r 6 or ^ 7 becomes doubtful. 

It is true that these higher orders in ^ do not need to be known very 
accurately, for they contribute but little to the final result ; nevertheless, an 
improvement in the calculation is called for. 

An easy escape from the difficulty is provided by the series for yfr n given 
fn (IV). When 7] is moderate, as here, and n somewhat large, the conver- 
gence is good, and we may calculate ty e and ty% by a perfectly straight- 
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forward process to a high degree of accuracy. Having obtained ^ 6 and ^ 5 , 
we may then use the sequence formula (YI) in the reverse direction to 
deduce in turn ^ 4 , \Jr 3 , etc., without loss of accuracy. The values of ^ n and 
i|r n , calculated for the various cases, are here tabulated. It is possible that 
they may prove useful in kindred enquiries. From them are ' deduced the 
complex values of M and N as already explained. In this case the tables 
give the logarithms of the real and imaginary parts. 







Tables of ¥ 


and ty. 




n. 


V. 


*. 


V. 


$, 


7) = 


* 1. 


7] = 


1-5. 



1 
2 
3 
4 
5 
6 


-54CT30 
4 -1453 
54 -075 
1747 '5 


-84147 
-90351 
-93053 
-94569 
-95541 


l 

-04716 
1 1 -3929 
8 -9713 
117 -89 
3049 -8 


-66499 
-79235 
-84900 
-88121 
-90204 
0-91664 
-92743 




v =, 


2 -25. 


7) = 


27/8. 



1 
2 
3 
4 
5 
6 


-0-279188 
+ -295636 
+ 1 -70318 
+ 9-73119 
+ 99 -9041 
+ 1933 -51 
+ 51793-7 


0*345811 
-577175 
-685527 
-749095 
-791075 
-820928 
-843272 


-0-288262 
-0-136838 
+ -199406 
+ 1-03318 
+ 4-61147 
+ 31 -1003 
| + 332 -545 

! 


-0-068531 
+ -238183 
+ -403904 
+ 0-509211 
+ 0-582438 f 
+ -636441 [ 
+ -677969 


7} - 


-2. \ 


7} = 


= 3. 



1 
2 
3 

4 
5 
6 


-0-208073 
+ 0-525919 
+ 2*75247 
+ 19 -4823 
+ 263 -495 
+ 6039 -32 
+ 206485 


-454649 
-653096 
-744180 
-796977 
-831564 
-856018 
-874244 


-0-329997 
-0-062959 
+ -445063 
+ 1 -97564 
+ 10 -7141 
+ 96*1227 
+ 1357 *05 


-047040 | 

-345677 

-497729 

-591312 

-655080 

'701448 ; 

-736733 




7) = 


7/4. 


V =5 


21/8. 



1 
2 
3 
4 
5 
6 


-0*101855 
+ -864129 
+ 4-73135 
+ 44-1968 
+ 811-861 
+ 24815*9 


-562278 
-725411 
-799021 
-841250 

-868721 
-888047 
-902393 


-0-331241 
+ -070827 
+ -875249 
+ 4 -08595 
+ 29 -3550 
+ 363 -049 


-188160 
-460481 
-592807 
-672133 
-725271 
-763446 
-792237 
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Tables of M and N. 



n. 


M„. 


JN«. 


1 

2 
3 


7] = 

[2 -4483]-* [4 -900] 
[5 -879] 
[5 -08] 


[1-2635]-* [2 -5424] 
[2 -0109] - i [I -0218] 

[4 -43] 


1 
2 
3 

1 4 


07 = 

' [1-85853]-* [2 -74171] 
[2-100] -*[4-200] 
[4 -633] 
L6 -99] 


3/2. 

[1*6160]-* [1-3389] 

[2 -8566] -i [3 -715] 
[3-622] -*[5-243] 

3 -16] 


1 
2 
3 

■ 4 
5 
6 


1 = 

[I -52056]-* T 1-94164] I 
[I -36224] - * [ 2 -749601 
[2 -18132]-* [ 3-36275] 
[I -87043]-* [ 7-7391] 
[5-377] -*[10-75] 
[7-80] -*[13-6] 


* 

9/4. 

[1 -69832]-* [ 1-72202] 
[I -58138]-* [ 1-24719] 
[2-77736J-*[ 3-55628] 
[3-6612] -*[ 5-322] 
[4-3277] -*[ 8-66] 
[6-88] -*[H-77] 


* 
■ 

' 1 
2 
3 
4 
5 
6 


1 = 

[I -68350]-* [1-80011] 
[2 -98851]-* [3 -98120] 
[3 -73029]-* [5 -46059] 

[4 -32695]-* [8 -654] 
[8 -7764] 
[7 -0937] 


[I -69493]-* [1-63560] 

[1-40882]-* [2 -84950] 
[2 -44577]-* [i -89188] 
[3 -23189]-* [6 -4638] 
[5 *8390] 
[8 -2901] 


i 

1 

2 
3 
4 
5 


V = 

[1-63910]-* [1 -40580] 

[2 -57067]-* [3 -14194] 
[3 -22163]-*' [8 -44] 
[5 -708] 
[6 -042] 


7/4. 

[I -67342]-* [I -52299] 
[1-16474]-* [2 -33907] 
[2 -06571]-* [1-1315] 

. [4 -7389] 
[5 -228] 



For a systematic calculation it appears preferable to employ the forms 
(k), (I) rather than the expansions in powers of ft. The functions involving 
P, once calculated for the various values of n, are then available for any 
value of 7). The following table exhibits the algebraic forms : — 
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n. 



1 
2 
3 
4 
5 
6 



(2 n + 1) P'„ 

w (n + 1) 



tA 4 

_7 
8 
9 



(-1 + 5 M 2 ) 
K-3^ + 7^) 
^(1-14 a* 8 + 21 /a 4 ) 
11(5^-30^ + 33^) 



(2n + l),P^ (2j> + 1)Pwt 
n (n + 1) 



1(1 



■ 2 ) 



(llju-15ju 3 ) 
■|(-3 + 27jk 8 -28/i 4 ) 
iJ(-29ju + 126/* 3 -105 J u 5 ^ 

(5 - 100 m 2 + 285 ^ - 198 /u 6 ) 



16 



The particular cases /a = 0, /a = ± 1 have already been considered. It 
will be observed that both functions are either entirely odd or entirely 
even. 

The tables give the logarithmic values for jj, = 0, J, |, f , 1. Thus, in the 
first table for yu, = J, ?i = 4, we have the number whose logarithm is T85776, 
and this is to be taken negatively when /j. = -f f. Moreover, since this 
function of /jl is odd, the sign is to be reversed, i.e., made positive, when 



/* = — i- 



Table for ( 2 ^ + 1 ) p/ » . 



w. 


fi = 0. 


u — i 
A* — 4' 


A* — i- 


fJL = f. 


fi = 1. 


1 1 

2 
3 

! 4 
! 5 

1 6 


+ [0-17609] 


- [I -94201] 



+ [I -83727] 



+ [0 -17609] 

+ [I -79588] 
-[1-77928] 
-[1-85776] 
+ [I -15331] 
+ [1 -82016] 


+ [0 -17609] 
+ [0 -09691] 
+ [1-33995] 

- [I -84703] 
-[1-91191] 

- [I -24977] 


+ [0 -17609] 
+ [0 -27300] 
+ [0-20029] 
+ [1-89819] 
- [I -19988] 
-[1-94131] 


+ [0 -17609] 
+ [0 -39794] 
+ [0 -54407] 
+ [0 -65321] 
+ [0 -74036] 
+ [0 -81291] 



Table for mP'»(^ + 1)_y 2 n + 1) p„. 



11. 


At = 0. 


A* ~ i- 


M = i- 


fi = f. 


. A* = 1. 


1 

2 
3 

4 
5 
6 




+ [0 -39794] 



- [0 -52827] 


+ [0 -60879] 


-[1-57403] 
+ [0 -33995] 
+ [0 -34265] 
r- [0 -20401] 

- [0 -56836] 

- [1 -17713] 


-[1-87506] 
+ [0-09691] 
+ [0 -50132] 
+ [0 -35218] 
- [0 -14504] 
- [0 -63256] 


- [0 -05115] 

- [1 -49485] 
+ [0 -22573] 
+ [0 -57335] 
+ [0 -64946] 

+ [0 -47638] 


- [0 -17609] 

- [0 -39794] 

- [0 -54407] 

- [0 -65321] 

- [0 -74036] 

- [0 -81291] 



The way is now clear for the last step, viz., the calculation of (yZ~-zY)fr 
and Y by summation with respect to n of the values given in (/) and (h). 
As to this there is little to call for special remark. Of course the odd and 
even terms are added separately. 
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Eigl 1. — Curves showing the moduli 
of Y (dotted) and (yZ-zY)jr (full) 
for various values of v\. The ab- 
scissa represents jx. 




+i-o 
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The tables exhibit the complex values and the modulus, subject to multi- 
plication by the factor (kr)~ l e iki ~ ct ~ T \ Thus for 77 = 1, //,=-— 1, we have 

-9/7 ?Y p ik(ct-r) p ik(ct~r-p) 

y/j ^ - e (0-3457-ix 0-0538)= 0'3499 e - — 



Icr 



Jcr 



where p influences the phase only. 

As regards the values of % or 7cE, or 2ir'R/X, the numbers 1, 3/2, 9/4 
were at first chosen. Subsequently it appeared desirable to fill up the gap 
at the more interesting stage, and so calculations for ?? = If , rj = 2 were added. 
The values of the modulus are plotted in the curves, fig. 1, as functions of pu 
for the various values of 97. It will be remembered that in all cases the 
refractive index is supposed to be 1*5. 



fX. 


yZ-zY 

• 


Modulus. 


1 1 

1 

i X . 

i 


Modulus. 


-1 


0*3457-ix0*0538 


1 = 
-3499 


= 1. 

0-3457-1x0*0538 


*3499 f 


3 

4: 


0*2525-ix 0-0404 


*2557 


0*3266-1x0*0534 


-3310 i 


1 

— 2 


-1669-1 x 0*0272 


*1691 


0*3083-ix0*0530 


-3128 i 


1 

4 


0'0884-ix0*0141 


0*0895 


0'2904-ixO -0527 


-2951 ; 





0-0165-ix0*0009 


*0165 


0*2731-1x0*0523 


*2781 


+ i 


~0*0490 + ix0-0121 


*0505 


0*2566-ix0*0519 


*2618 1 


+ ^ 


-0-1083 + ix0-0250 


*1111 


0*2405-ix0-0516 


-2460 


+ 4: 


-0-1619 + 1x0-0380 


-1663 


0-2250-ix 0-0512 


*2307 ! 


+ 1 


-0-2103 + 1x0*0508 


-2164 


0-2103~ix0*0508 

I 


*2164 

. .- — — — 1 . — . 






1 = 


1-5. 




-1 


1-1902-1x0-4237 


1 -2634 


I*1902-ix0*4237 


1 *2634 


'A 

4: 


0-8463-ix0-3302 


-9084 


1-0212-1x0-3992 


1 '0965 


L 

™ 2" 


0-5647-ix0-2401 


-6136 


0*8643-ix0*3751 


*9422 


i 

4 


0-3388-ix0-1533 


-3719 


0-7190-ix0-3510 


*8001 





0-1629-ix0-0698 


-1772 


0*5844-ix0*3270 


*6697 


+ 4 


0-0316 + ix0-0105 


-0333 


0*4602-ix0*3032 


'Ouii 


4- i 
+ 2 


-0-0597 + ix0-0875 


-1059 


0*3453-ix0*2793 


*4441 


4. s 

T^ 4 


-0-1159 + ix0*1614 


-1987 


0*2390-1x0*2556 


*3499 


+ 1 


-0*1414 + 1x0-2321 


*2718 


0*1414-ix0*2321 


*2718 






V = 


1*75. 




-1 


1*8683-1x0*9404 


2-0917 


l*8683-ix0*9404 


2 *0917 


3 
4 


1-2802-1x0*7661 


1 *4919 


l*4992-ix0*8280 


1 -7127 


1 
2 


0*8330-1x0-6059 


1 -0301 


1-1668-1x0*7165 


1 *3693 


1 

. 4 


0-5081-ixO-4597 


'6852 


0*8694-1x0*6060 


1*0598 





0-2885-ix0-3273 


*4363 


0-6041-ix0'4964 


'7819 


+ i 


0-1592-ix0-2085 


*2623 


0*3682-1x0*3878 


*5347 


+ * 


0-1067-ix0-1032 


*1484 


0*1594-ix0 -2801 


*3223 


+ 1- 


0-1191-ix0-0112 


*1196 


-0 '0251-?; x 0-1734 


*1752 


+ 1 


0-1858 + ix0-0675 


-1977 


-0*1858-1x0-0675 


*1977 

„ — , — ™„ __ — —————— '■■■—■' — > 
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fj. 


#Z-zY 

• 


Modulus. 

i 


Y. 


Modulus. ,j 


-1 


2-4767-ixl-7984 


v ■- 

3 -0608 


= 2. 

2*4767-^x1-7984 


3 -0608 


3 

4 


l*4988-ixl-47l6 


2 -1005 


1-8327-^x1-4949 


2 -3650 


—.1 

a 


0*8053-ixl -1919 


1 -4384 


l-2909-txl-1980 


1 *76U 


— 4 


0*3478-ix0'9583 


1 -0195 


0-8417-^x0-9076 


1 -2379 





-0839-1 xO '7699 


-7745 


0'4759-ix 0*6235 


-7844 


+i 


-0 -0228-1 x 0*6257 


-6261 


0*1855-^x0 -3459 


-3925 


+4 


-0 '0040-1x0 -5247 


'5247 


-0-0372-^x0-0747 


-0834 


+ f 


+ 0*1132-1x0 -4660 


'4795 


-0*1988 + ix 0*1901 


-2751 


+ 1 


+ 0*3056-1x0 -4487 


-5429 j 


-0-3056 + ^x0-4487 


-5429 






V = 


2*25. 




-1 


3-0635-ix2*6979 


4 '0822 


! 3*0635-^x2-6979 


4-0822 


'A 

4 


1-5157-1x2 -0593 


2 -5569 


1 -9791-ix2-1286 


2-9065 


2 


-4857-ixl-5448 


1 -6193 


1-1464-^x1-5972 


1*9661 


-i 


-0-1392-^x1-1497 


1-1581 


0-5364-^x1-1036 


1 -2270 





-0 '4540-ix0-8696 


-9810 


0-1234-ix0*6473 


-6590 


+ i 


-0'5383-ix 0*6999 


'8830 


-O'1157-ix 0*2281 


-2558 


+ £ 


-0-4585-^x0-6364 


-7844 


-0-2018 + ixO -1544 


-2541 


+ * 


-0-2689-^x0-6745 


-7261 


-0-1531 + ix0-5003 


-5232 


+ i 


-0 -0139-?; x 0-8101 


-8102 


+ 0-0139 + 1x0-8101 


0-8102 



When 7) = 1, i.e., when the circumference of the sphere equals the wave- 
length in the general medium, there is comparatively little departure from the 
conditions appropriate to the infinitely small sphere. When the secondary 
ray is in a direction perpendicular to that of primary vibration, viz., when 
y = 0, the vibration (Y) varies but little with /n. It is otherwise when the 
secondary ray is in the plane containing the directions of the primary ray 
and of primary vibration. The second column shows that in that case the 
secondary vibration varies greatly with \x, and that in the perpendicular 
direction (/& == 0) the intensity is very small. This is the direction in which 
there would be no effect, however large the sphere, if the square of the 
•difference of optical quality could be neglected absolutely. 

The departure from the law of the infinitely small sphere becomes more 
apparent when ?/ = 1*5. The principal feature here is the displacement of 
the maximum polarisation from fi = to the neighbourhood of fju = + \, i.e., 
to a direction inclining hack/wards along the course of the incident light. 
When the incident light is unpolarised, the diffracted light is polarised at all 
angles (other than /u = + 1), and in the same sense as from infinitely small 
particles. 

The latter law ceases to hold when rj rises to about 1*75. A third neutral 
point (direction of no polarisation) enters at jjl = + 1 and moves inwards. 
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When rj = 2, the change has proceeded so far that the neutral point is in the 
neighbourhood of /* = 0, the polarisation for positive fi being the reverse of 
that, found for infinitely small particles. At fi = + -£ the reversed polarisation 
is very pronounced. When r\ attains 2*25, the neutral point is still further 
on (near //, = — \), and the maximum reverse polarisation is near //, = ^. It 
will be seen that these changes occur very rapidly in the neighbourhood of 
rj = 2. 

It would be possible to follow these calculations to greater values of rj, and 
such an extension would not be without interest, but the arithmetical work 
would soon become heavy. Also, without increasing ??, the refractive index 
might be varied. Some data for such calculations are included in the tables 
already given. 

So far the primary light has been supposed to be either unpolarised, or, 
if polarised, to be plane polarised in one or other principal plane. In the 
more general case recourse must be had to the complex expressions, which 
allow the effects of the two principal components of the primary light to be 
combined, with the necessary allowance for any phase-difference. But it is 
scarcely necessary to pursue the treatment of this part of the subject. 

Experimental. 

So far as I know, the only experimental work bearing at all closely upon 
the preceding calculations is that of Govi and Tyndall (' Proc. Roy. Soc./ 
January, 1869), who observed the light scattered from smoke of various kinds. 
Tyndall established the existence of neutral points at angles dependent upon 
the condition of the smoke, but his indications are not very precise. For 
example, he does not distinguish between density due to increased number 
of particles in a given space and that due to increasing size, of which the 
effects may be expected to be very different. Nor, again, do I see any 
suggestion that the position of the neutral point varies with the colour of 
the light. Nevertheless, TyndaU's description is suggestive, and it formed 
the starting-point of my own work upon the subject. 

In continuation of the experiments already referred to at the commence- 
ment of this paper, I have lately repeated the observation on the nascent 
precipitate formed when a little acid is added to a dilute and well-filtered 
solution of hyposulphite. The liquid was contained in a small beaker placed 
in a dark room and exposed to a beam of sunlight. The observation of the 
light dispersed at various angles may be made with a nicol, but a double 
image prism is better. This may be mounted at one end of a tube of which 
the other carries a small aperture of such dimensions that two oppositely 
polarised rectangular fields are seen in close juxtaposition. In the early 
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stage the double-image prism, suitably rotated, shows complete polarisation of 
the light scattered along BO, fig. 2, at right angles to the primary rays AB, 
Somewhat later the minimum light seen along BC shows TyndalFs " residual 




Fig. 2. 



blue/' indicating that for this colour the polarisation is no longer complete. 
But if the direction of observation be changed towards BD, the polarisation 
is improved. As the particles grow, a stage is reached in which the effects 
for blue and red light are well contrasted. If a blue glass be inserted in the 
course either of the primary or of the secondary light, that polarised com- 
ponent which originally vanished (along BC) is found along BD to be the 
brighter of the two, showing that the polarisation is reversed The substi- 
tution of a red for the blue glass restores the original character of the 
polarisation, the increase in the wave-length compensating the growth of the 
particles. It is not pretended that the observed contrast is as great as the 
theoretical curves might suggest ; nor could this reasonably be expected, seeing 
that the particles are not all of the same size. 

April 30. — I have recently learned from Prof. Love that he had become 
aware of the error in equation XVIII, and that a correction was published in 
Math. Soc. Proc./ vol. 31, p. 489. 



